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o , 

^ O ^ The description of static critical phenomena in terms of finite size scaling (FSS) relations, developed by Fisher et al. 

[Q,^ has been extended by Halperin, Hohenberg, Ma and Suzuki to include dynamical properties of the system. 
^ Later, Janssen et al. and independently Huse [|i found evidence for an universal behavior far from equilibrium. 

00 ; 

00 ' As discussed by Janssen et al. one finds an universal behavior already in the early stages of the relaxation process for 
,, systems prepared at an initial state characterized by non-equilibrium values of the order parameter. As a consequence, 

' I I they could advance the existence of a new critical exponent 9, independent of the known set of static exponents and 
' of the dynamic critical exponent z. This new exponent characterizes the so called "critical initial slip", the anomalous 
d increasing of the magnetization when the system is quenched to the critical temperature Tc. 

^ , That new universal stage has been exhaustively investigated to confirm theoretical predictions and to enlarge our 
^ knowledge on phase transitions and critical phenomena. In this sense, several models and algorithms [Q-^ have 
O ■ been used, as toy models, in order to check the ability of the new approach in obtaining dynamic and static critical 
^ ' exponents. Results are in good agreement with pertinent results for static exponents and seems to be confident even 

for the new critical exponent d. However, a reliable technique to obtain the dynamic exponent z is lacking. A first 
^ . proposal by Li et al. |l0| using a time-dependent Binder's cumulant yields estimates with low precision when compared 

with other techniques ||ll|,|l^. An alternative way which uses another kind of cumulant, proposed by Zheng |13|, gives 

the right answer for the 2D Ising model but fails in determining the value of z for the 3-state Potts model [|l3| and 

for the Ising model with multispin interactions [p^ . 

In this letter, we introduce and check a new technique to obtain the exponent z, combining the behavior of the 
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order parameter and its second moment when the system is submitted to different initial conditions. 
Before presenting our proposal, we shall review the main results in short time dynamics. 

Although Halperin et al. |^] have studied systems with different dynamics, we consider only systems without 
conservation laws, the so called Model A because our discussion aims dynamics generated by heat bath dynamics. 
Therefore, we consider a magnetic system prepared at high temperature (T >> Tc) with a small nonzero magnetization 
mo, (this can be achieved with a small external magnetic field h) and quenched to the critical temperature Tc without 
any external magnetic field. If the system is allowed to relax towards equilibrium with the dynamics of model A, the 
magnetization obeys the following scaling relation (generalized to fcth moment). 

Here h is an arbitrary spatial scaling factor, t is the time evolution and r is the reduced temperature, r — {T ~Tc)/Tc. 
The exponents /? and v are the static critical exponents, while z is the dynamic one. are the fcth 

moments of magnetization. This scaling relation depends on the initial magnetization mo and gives origin to a new, 
independent critical exponent a;o; the scaling dimension of the initial magnetization, which is related to 6. 

From Eq. (^ we can derive the power law increasing of the magnetization, observed in the initial stage of the 
dynamic relaxation. For this purpose, we consider large lattice sizes L at r = with b = t^^^. This leads to the 
scaling relation 

M{t, mo) = t-'^'''m{l, t^°/^mo) (2) 

for the first moment of the magnetization M^^^ = M . By expanding this equation for small mo, we have the following 
power law, 

M{t) - mo*^ , (3) 

where, as anticipated, we identify 6 = (xo — P/i')/z. Here we also have the condition that t^°^^mo is small, which 
sets a time scale to ~ mg where that phenomena can be observed. 

On the other hand, it has been realized the existence of another important dynamic process from an initial ordered 
state p8|-p0| , which represents another fixed point in the context of renormalization group approach. This leads to a 
different universal behavior of the dynamic relaxation process also described by Eq. (^ with mo = 1. In particular, 
dealing with large enough lattice sizes at the critical temperature (r = 0), one obtains a power law decay of the 
magnetization 

M{t) - t-f^/"' , (4) 

which follows from Eq. (^ when we choose b^'t — 1 and shows the average magnetization is not zero as would be 
expected from a disordered initial state. 

Equation (|l|) and their particular forms in (|^) and (^ can be used to determine relations involving static critical 
exponents and the dynamic exponent z p"3|Jl^ ]. 

The observables in short-time analysis are described by different scaling relations according to the initial magneti- 
zation. In particular, the second moment M^'^\t^L), 



(aia-i) , (5) 

behaves as L"*^ since in the short-time evolution with initial condition mg ~ the spatial correlation length is very 
small compared with the lattice size L. Thus we have 

M(2) {t, L) = t-^^'''' Af(2) (1, i-i/U) - t^d-^P/-)/^ . (6) 

Moreover, under this initial condition one can also define the time-dependent Binder's cumulant at the critical tem- 
perature, 

which leads dXT — Tc to the FSS relation 

U{t,L) ^U{b-'t,b-^L), (8) 

and the exponent z can be independently evaluated through scaling collapses for different lattice sizes [p^pT| . 

In order to obtain more precise estimates for the dynamic exponent z, another cumulant has been proposed 
It is given by 

and should behave as 

U2it)^t^/' (10) 

when one starts from an ordered state. In this case, curves for all the lattices lay on the same straight line without 
any re-scaling in time and results in more precise estimates for z. However, application of this procedure has not been 
successful in at least two well known models: the two-dimensional q ~ 'i Potts model and the Ising model with 
three spin interactions in just one direction p^ . The reason for the above disagreement could be related to the value 
of the second term of r.h.s in Eq. when mo — 1. 

A plausible way to circumvent this problem is to work with different initial conditions. For this we decided to 
follow the evolution of the ratio F2 — Af '■^^ /M^ but using different initial conditions to calculate each one of the mean 
values. The reason is we know the behavior of the second moment of the magnetization when samples are initially 
disordered (mo — 0) and also the dependence on time of the magnetization of samples initially ordered (mo = 1). 
Under the above mentioned conditions the ratio behaves as 

which has the same potential law that the cumulant mentioned before but requires two independent simulations 
instead of one used for calculating U2- 

In short-time Monte Carlo (MC) simulations, the time scale t is settled in units of whole lattice updates, and does 
not depend on the initial conditions. However, the same dynamics (Metropolis, Glauber or heat-bath) should be used 
in both simulations @J^. 



Now we present our estimates obtained for 2D Ising model, q = 3 and q = A Potts models. 

We have performed independent heat-bath (HB) MC simulations for a large lattice (L = 144) according to the 
required initial conditions and the evolutions have been done until the maximum time t — 200 MC sweeps with 
N = 10000 samples. This gives averages for the magnetization and its second moment. We have performed this kind 
of simulation 20 times to obtain our final estimates for each moment in function of t. Since all runs are independent 
due to the choice of random numbers, we have a total of 400 time series for Eq. (^ij). 

Our results are shown on a log-log scale for the time interval [10, 100] in Fig.l and Fig. 2, respectively for three and 
four-state Potts model (a similar figure is obtained for 2D Ising model) and they appear to be clearly consistent with 
a power law in that time interval. The error bars are also presented in those figures but they are hardly seen in that 
scale. 

Our estimates for z comes from a least square fitting in the time interval [ti, tf]. Due to our small statistical errors, 
we can make a systematic study for the range in t where we find acceptable goodness-of-fit Q |^l[. As examples, we 
obtain for 2D Ising model, z = 2.1435(2) in the time interval [10, 200] with Q = 10"250, z = 2.1359(3) in [50, 200], with 
Q — 10^^^, and the most acceptable Q — 0.99 in [30, 90], which yields z — 2.1565(7). This value is presented in Table 
1, where we also include some estimates for comparison. Here we observe that our result is in agreement with more 
recent estimates within two standard-deviations. This indicates our statistical errors are presumably underestimated 
possibly due to corrections to scaling. 

We complete the overview in Table 1 with data from [^. Estimates have been obtained from the long time 
behavior of the magnetization for the square lattice (sq), z = 2.168(5), for the triangle (TP), z — 2.180(9), and for 
the honeycomb (he) lattice, z ~ 2.167(8), while from damage spreading in short-time the authors quote 2.166(7), 
2.164(7) and 2.170(10) for sq, TP and he lattices. 

For q — 3 Potts model, our study monitoring Q gives z = 2.198(2) in the interval [50, 90] with Q — 0.82. This 
result agrees with z = 2.203(11) (Ref. jl^]) obtained from the Binder cumulant in Eq. (0), and with the estimate 
obtained in j|] with HB algorithm from the second moment M^^^t) in short-time analysis. 

The value in [^ refers to TP lattice, presenting further numerical evidence (comparing [^ and ^) to the dynamic 
universality. Reference [jl^ also presents the value 2.14(3), obtained from Eq. (^), in clear disagreement with 2.203(11) 
as commented in [p^ . On the other hand, our estimate from Eq. ([ll| ) gives a value in full agreement with the Binder 
cumulant and collapse data analysis in short time. 

The case q — ^ has been less studied. Our analysis gives z ~ 2.290(3) in the interval [60, 90] with Q = 0.72. Here, 
we stress the importance of monitoring Q, since we may find values for z as large as z = 2.3483(2) in [10, 200] but 
with unacceptable value, Q = lO^^eg pj. z = 2.3532(3) (2) in [10, 100] with Q = lO-^o^. 

As a final comment, the suspected z as been weakly [ p8| or even independent on q [ pTj is not supported by the most 
recently results presented in Table 1. 
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Figure Captions: 



Figure 1. Time evolution of F2{t) for 2D three-state Potts model. 
Figure 2. Time evolution of F2{t) for 2D four-state Potts model. 



TABLE 1. Dynamic exponent z for 2D Ising model, three and four-state Potts model. 



Reference (year) 



Ising 



This work 
(2000) 
(1999) 
(1998) 
(1998) 
(1998) 
(1997) 
(1997) 
(1995) 
(1995) 
(1993) 
(1992) 
(1987) 
(1986) 



2.1565(7) 

2.1667(5) 

2.153(2) 

2.153(4) 

2.16(2) 

2.137(8) 

2.166(7) 

2.155(3) 

2.172(6) 

2.16(4) 

2.165(10) 

2.16(2) 

2.16(5) 



2.198(2) 

2.191(6) 

2.203(11) 

2.14(3) 



2.196(8) 



2.16(4) 
2.43(15) 



2.290(3) 



2.18(3) 
2.36(20) 



1 2D TP lattice. 

^ from scaling collapse (Table 3). 
3 applying Eq.(|). 

^ from HB algorithm, while z = 2.137(11) andz = 2.198(13), respectively for Ising and g = 3 model from Metropolis 
algorithm. 




(Fig.l) 1. Time evolution of F2{t) for 2D three-state Potts model. 




(Fig.2) 2. Time evolution of F2{t) for 2D four-state Potts model. 



